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Abstract 



In this thesis, we have investigated the possibihty of large cross sections at large 
multiplicity in weakly coupled three dimensional 0^ theory using Monte Carlo Sim- 
ulation methods. We have built a Beowulf Supercomputer for this purpose. We use 
spectral function sum rules to derive a bound on the total cross section where the 
quantity determining the bound can be measured by Monte Carlo simulation in Eu- 
clidean space. We determine the critical threshold energy for large high multiplicity 
cross section according to the analysis of M.B. Volosion and E.N. Argyres, R.M.P. 
Kleiss, and C.G. Papadopoulos. We compare the simulation results with the per- 
turbation results and see no evidence for large cross section in the range where tree 
diagram estimates suggest they should exist. 
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Overview 



Chapter 1 of this thesis is the introduction to the high multiphcity cross section 
problem and the idea to solve the problem using Monte Carlo simulation. Chapter 2 
describes the detail of design and construction of the Beowulf parallel supercomputer. 
In chapter 3, the sum rules and the bound for the cross section are derived. 

In chapter 4, a perturbative estimate of the critical energy is given. And the 
relationship between the critical energy and the Euclidean momentum on the lattice 
is determined. In chapter 5, the Monte Carlo simulation is described. The scaling 
behavior and continuum limit are also discussed. In chapter 6, I will use the pertur- 
bation method to calculate the cross sections. Chapter 7 presents the results from 
the simulation. Finally, I will compare the results with the perturbative calculations 
and discuss conclusions. 
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Chapter 1 
Introduction 



In 1990, Ringwald and Espinosa presented a calculation of possible instanton 
induced baryon number violation accompanied by enhanced production of large num- 
ber of Higgs and W,Z bosons at high energy(£' ~ 72/(7, where ^ is the number of out 
going particles and g is a. generic coupling constant). This has called attention to a 
more general and profound problem of quantum field theory. The problem is: for a 
given renormalized quantum field theory with weak renormalized coupling constant, 
what is the nature of the multi-particle production in that theory. One can think two 
extreme cases. In the first case, the cross section is perturbatively small[~ x 0(1)]. 
In the second case, when the number of out going particles n ~ 1/(7, the multipar- 
ticle production cross section becomes unsuppressed and the inclusive cross sections 
saturate unitarity bounds. 

It has long been known that the perturbation-theory expansion in a theory with 
weak coupling fails in high orders because of the factorial growth of the coefficients 
in the series 0. This problem can not be solved perturbatively because the pertur- 
bation series diverges rapidly at high order. In fact, independent of the instanton 
calculations, Goldberg |^, Cornwall ||^ and, later, Voloshin^, Argyres, Kleiss, and 
Papadopoulos point out that in 0'* scalar field theory the contributions of just the 
tree graphs to the multiparticle production amplitudes gives the nig"' behavior. This 
suggests the possibility of the second case. 

Here we have a different approach to solve this nonperturbative problem. We pro- 
pose to use Monte Carlo simulation on a lattice to determine an upper bound on the 
cross section of the multiparticle production in 01 theory. The original application 
of this approach was to two dimensional 0| theory [^]. However, the case of 1 + 1 di- 
mensions is a very special case in field theory. (no angular momentum... etc) Although 
the invesgation in 1 + 1 dimension showed no signal of large multiparticle production, 
the result might be special for the 1-1-1 dimension case. In this research, we will 
use Monte Carlo simulation to establish a bound on the multiparticle cross section in 
three dimensions(l -|- 2 dimensions in Minkowski). The tree diagram considerations 
which suggest large 1 ^ n amplitudes in 0^ are mainly combinatoric, and hold as 
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well in three dimensions as in two or four. Working in three dimensions has several 
advantages. First, we can work on large A^'^ lattices as long as our computers can pro- 
vide the memory and give us reasonable performance. Second, three dimensional 0"^ 
theory is a superrenormalizable theory. The wave function renormalization constant 
Z is finite in the 0| theory. The other important issue in our simulation is that we can 
do the simulation in three dimensional Euclidean space. The constructive quantum 
field theorists have proved that the (^\ field theory exists as the analytic continuation 
of the continuum limit of the Euclidean lattice theory |jlO|, [ill. In four dimensions, 
nonperturbatively, we would have to deal with the presumed triviality of the theory. 

Even though we worked on three dimensional Euclidean lattices rather than four, 
the simulation itself is still a huge computing task. This is because the many degree 
of freedoms required for sensitivity to multiparticle production requires large lattice. 
For this reason, we will still need a supercomputer to do the job. To solve this 
problem, we built a Beowulf parallel supercomputer with sixteen CPUs. We also 
use the Message Passing Interface(MPI) function call to parallelize our simulation 
program. We are able to achieve overall 7.4 Gflops at High Performance Computing 
Linpack Benchmark(HPL). The theoretical peak performance for sixteen CPUs is 12 
Gflops. Compare with 16 Gflops CRAY C90 at 16 CPUs, our parallel cluster has a 
much better price/performance ratio. 



Chapter 2 

Beowulf Parallel Supercomputer 



2.1 The Beowulf Project History 

The Beowulf Project was started by Donald Becker when he moved to CESDIS in 
early 1994. CESDIS was located at NASA's Goddard Space Flight Center, and was 
operated for NASA by USRA|12[. 

In the summer of 1994 the first Beowulf 16 node demonstration cluster was con- 
structed for the Earth and Space Sciences project, (ESS). The project quickly spread 
to other NASA sites, other R&D labs and to universities around the world. The 
project's scope and the number of Beowulf installations have grown over the years. 
Today, the fastest Beowulf cluster can reaches 143 Gflops (HPL benchmark) with 
528 Pentium III 800 CPUs which ranks at 126 in the topSOO supercomputer list in 



Nov, 2000|]13|. With more and faster CPUs and even faster network links, Beowulf 
supercomputer can compete with any commercial high performance supercomputer 
such as SGI Origin 2000 or Cray T3E and costs much less. 



2.2 Definition of the Beowulf Supercomputer 

Beowulf is a multi computer architecture which can be used for parallel computations. 
It is a system which usually consists of one server node, and one or more client nodes 
connected together via Ethernet or some other network. It is a system built using 
commodity hardware components, like any PC capable of running Linux, standard 
Ethernet adapters, and switches. It does not contain any custom hardware compo- 
nents and is trivially reproducible. Beowulf also uses commodity software like the 
Linux or FreeBSD operating system. Parallel Virtual Machine (PVM) and Message 
Passing Interface (MPI). The server node controls the whole cluster and serves files 
to the client nodes. It might also be the cluster's console and gateway to the outside 
world if there is no particular annex machine. Large Beowulf machines might have 
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more than one server node, and possibly other nodes dedicated to particular tasks, for 
example consoles or monitoring stations. The client nodes in the cluster don't sit on 
people's desks; they are dedicated to running cluster jobs only. In most cases client 
nodes in a Beowulf system are dumb and are configured and controlled by the server 
node, and do only what they are told to do. In a disk-less client configuration, client 
nodes don't even know their IP address or name until the server tells them what 
it is. One of the main differences between Beowulf and a Cluster of Workstations 
(COW) is the fact that Beowulf behaves more like a single machine rather than many 
workstations. In most cases client nodes do not have keyboards or monitors, and are 
accessed only via remote login or possibly serial terminal. Beowulf client nodes can 
be thought of as a CPU + memory package which can be plugged in to the cluster, 
just like a CPU or memory module can be plugged into a motherboard. In general, 
the user can only login to the front end and submit their jobs to the cluster. The 



queueing system will receive the jobs and run them in different nodes|14 



To get better performance, one can also use some dedicated parts which can not be 
purchased on the commodity market. This is true especially when people need a very 
fast link between the nodes. One can use Giga bits Ethernet or Myrinet(240MB/s at 
64-bit, 33MHz PCI) or even Quadric(340MB/s) system for the fast link. The other 
way to improve the floating point performance is to use DSP chips such as Texas 
Instrument TMS320C6xxx. These chips can boost the floating point performance up 
to 3 Gfiops. However, the drawback is that they are expensive and some of them lack 
driver support. 



2.3 Architecture and System Design 

Our Beowulf Supercomputer is designed to solve large-scale scientific problems. Be- 
cause we need to run our simulation on a large lattices, we also need huge amount of 
memory. Each node within the cluster is running the freely available Linux operat- 
ing system. Parallel capabilities are achieved through the use of the MPI and PVM 
communication protocols. Job scheduling and execution for all users is handled by 
the Distributed Queueing System(DQS). The system contain a server node, a annex 
node, and 16 client nodes. All of these machines are networked by 100 BaseT Ether- 
net with low latency switch, (see Figure p.l| ) This architecture is so call distributed 
memory(or local memory) architecture which is good to run both Single Program, 
Multiple Data(SPMD) and Multiple Instruction, Multiple Data(MIMD) program. 
The machine was assembled entirely from individually chosen commodity hardware 
components. With the exception of the Ethernet switch, all of the computer hard- 
ware within the cluster was built from the component level. This low-level approach 
allowed the developers to optimize the overall performance by providing us the free- 
dom to choose the fastest and most reliable components to fit our budget. Likewise, 
all matters of OS and software package configuration were carried out in- house in an 



effort to achieve optimal settings for speed and stability 1 15 . 
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Beowulf Parallel Cluster 



Department 
Network 




node 1 node 2 



node 16 



Figure 2.1: The architecture of the Beowulf Supercomputer 
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2.3.1 Hardware Configuration 

The server node contains: 

• Pentium III 500 MHz CPU 

• ASUS P2B-S Motherboard 

• 384 MB Samsung SDRAM 

• 27 GB on two Ultra2 Wide LVD Quantum Atlas III SCSI HDs 

• Seagate 12/24 GB DAT Tape Backup System 

• Kingston KNE-IOOTX Tuhp-based Ethernet card 

• PCI video card 

• Linksys I/O console selector switch 
Each client nodes has: 

• Pentium III 744 MHz Coppermine CPU 

• ASUS P2B-F Motherboard 

• 768 MB Micron PC133 SDRAM 

• 4.3 GB UDMA Quantum Fireball HD 

• Kingston KNE-IOOTX Tuhp-based Ethernet card 

• PCI video card (for boot-up only) 

The annex machine is the gateway of the cluster. It has: 

• ABIT BP6 Dual Celeron Motherboard with Two 450MHz CPUs 

• 256 MB Samsung SDRAM 

• 6 GB UDMA33 IDE HD 

• Kingston KNE-IOOTX Tuhp-based Ethernet card 

• 3Com lObase-T Ethernet card (for outside LAN) 

• PCI video card 
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2.3.2 Software Configuration 

All of the nodes are running Linux Operating System. The Clients nodes share the 
user and utility files with the other via NFS (network file system). For the mes- 
sage passing protocol, we use message passing Interface (MPI), Parallel Virtual Ma- 
chine(PVM), High Performance Fortran(HPF), and newly developed Cactus Code. 
The DQS queueing system can handle both serial and parallel jobs. 

2.4 Performance and Benchmark 

The performance of our Beowulf Cluster is quite remarkable. We tested it with two 
different parallel benchmark program which are recommended by the topSOO cluster 
web site, NAS and HPL0. The NAS Parallel Benchmarks (NPB) are a set of 8 
programs designed to help evaluate the performance of parallel supercomputers. The 
benchmarks, which are derived from computational fluid dynamics (CFD) applica- 
tions, consist of five kernels and three pseudo-applications. 

The detail description for each benchmark routine is following: 

• CG - solution of a structured sparse linear system by the conjugate gradient 
method 

• FT - this benchmark contains the computational kernel of a three dimensional 
FFT-based spectral method 

• MG - Multigrid Benchmark uses a multigrid method to compute the solution 
of the three-dimensional scalar Poisson equation 

• LU - a simulated CFD application which uses symmetric successive over-relaxation 
(SSOR) to solve a block lower triangular-block upper triangular system of equa- 
tions resulting from an unfactored implicit finite-difference discretization of the 
Navier-Stokes equations in three dimensions 

• IS - parallel sort over small integers 

• EP - embarrassingly parallel benchmark 

• SP and BT - SP and BT are simulated CFD applications that solve systems 
of equations resulting from an approximately factored implicit finite-difference 
discretization of the Navier-Stokes equations. 

• BT code solves block-tridiagonal systems of 5x5 blocks; the SP code solves scalar 
pentadiagonal systems resulting from full diagonalization of the approximately 
factored scheme. 
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Table 2.1: NAS Benchmark 



Benchmark 


Number of 


Beowulf 


OilAi ioil/iZUU 


iN ame 


iNOQeb 




I'frnm M A Q M A Q A 
l^IIOm iN/\0,iN/\0/\J 


BT class=C 


16 


1250.88 Mfiops 


1088.3 Mfiops 


IS class=C 


16 


28.83 Mfiops 


26.9 Mfiops 


CG class=C 


16 


211.63 Mfiops 


148.7 Mfiops 


EP class=C 


16 


52.56 Mfiops 


54.8 Mfiops 


LU class=C 


16 


1380.13 Mfiops 


1169.5 Mfiops 


MG class=C 


16 


1071.59 Mfiops 


1616.8 Mfiops 


SP class=C 


16 


871.19 Mfiops 


800.4 Mfiops 


FT class=C 


16 


392.64 Mfiops 


NA 



Here we show the class C (large size) results in Table |2]1|. Compare with CRAY 
T3E 1200, one can see that the Beowulf Supercomputer is even faster then T3E 1200 
at same number of CPUs. However, for a heavy communication job, the T3E will be 
faster because of the faster interlink between the CPUs. 

The second benchmark is High Performance Linpack. This benchmark is the one 



used in the Top500 supercomputer list [ 13 1. HPL is a software package that solves a 



(random) dense linear system in double precision (64 bits) arithmetic on distributed- 
memory computers. It can thus be regarded as a portable as well as freely available 
implementation of the High Performance Computing Linpack Benchmark. Here we 
run three different sizes of the problem with different number of CPUs. The table 



shows the results. One can see the fastest result is 7.4 Gfiops at 16 CPUs(4 x 4). 



The Figure 0| shows that up to 16 CPUs, the speed scales up linearly. At the scale 



of 32 nodes or above [|T2|, the performance saturates the bandwidth of the lOOMbits 
network. The Giga bits network can boost the performance more. However, the 
Giga bits network suffers from the TCP/IP protocol which was designed for lOMbits 
network a longtime ago. For the Giga bits Ethernet link, one can only optimize either 
for the bandwidth or for the latency. 



Table 2.2: HPL Benchmark 



Matrix Size 


Block size 


P 


Q 


Time 


Results(Gfiops) 


20000 


200 


2 


4 


1501 sec 


3.55 


25000 


200 


2 


5 


2324 sec 


4.48 


34000 


200 


4 


4 


3545 sec 


7.4 
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HPL Benchmark Results 




Number of CPUs 

Figure 2.2: The HPL Benchmark shows the speed scales up Unearly. 



Chapter 3 

The Cross Section and Bound 



In this chapter, we want to determine the upper bound for the inclusive cross section 
of the multiparticle production. We will derive the bound in four steps in following 
four sections. First of all, we will specify the scattering process and define the inclusive 
cross section. Secondly, use the spectral representation of the two-point function and 
find the relation between the inclusive cross section and the spectral function. Thirdly, 
define G~^, Z', and m' and give the spectral representation of the function. In 
the last section, we will find the upper bound of the inclusive cross section. 



3.1 Definition of the inclusive cross section 

As we mentioned in chapter one, we want to study the amplitudes for one off- 
shell (time-like) (f) particle to go to n on-shell (p particles (see Figure |3.1|) . There are 
the amplitudes for which extensive tree diagram analysis have been done@, ||, H, |^, |]. 
These amplitudes may be considered as "physical" if we introduce a weak Yukawa 
coupling of the field to a fermion 

Cy = gyij^xp, (3.1) 
and compute the cross section for // to annihilate into n0 particles. (See Figure p.2|) 
This is similar to the e^e~ annihilation into hardrons in QCD. The transition 



probability for Figure |3.1| is 

,2 



m 



T{n,E) = - fd^n{E)\A{n;E)\^ (3.2) 
nl J 

where E = ^/s. The amplitude A{n] E) is 

A{n; E) = (pi, P2, ■ ■ ■ , Pn) (3.3) 

which is the connected (n+l) point Green function. The underscore means amputated 
Green function (no external line factors) and the slash means renormalized Green 
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n (]) 



Figure 3.1: The amplitude A{n\ E) for cf) —>■ ncf). 




function (no Z factions in the Lehmann-Symanzik-Zimmermann reduction formulas). 



It is convenient define a partial "cross section" of Figure 3.1 as a kinematic factor 



times equation ^ 



m 
2tts' 



;T{n;E) 



(3.4) 



for n out going particles. The the cross section of // — n0 in three dimensions 
then can be written as(to lowest order in Yukawa coupling but to all orders in the 
self- coupling]^) 



cTj/^„^(s) = -^\s(^m\ X a„(s) 



(3.5) 



Finally, we include all possible number of out going particles and define the inclusive 
cross section as 



a{s) 



II ^n( 



(3.6) 
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j^E;^/^^n|^(n+l)(Pl,---,Pn)r (3-7) 



27rs „ 



3.2 Spectral Representation and Sum Rules 

In the last section, we defined the inclusive cross section for the multiparticle final 
state. One can relate the inclusive cross section to the spectral function of the two 
point function. Let (poi^) be the canonical, unrenormalized scalar field in Minkowski 
space. We have the Wightman two point function 

W{x,y)^{n\cf>{x)cj>{y)\n) (3.8) 

The two point Green function (or r function) is the time ordered Wightman function 

Gix,y) = {n\TicP{x)cj>iy)\n) (3.9) 

Define the spectrum of physical state in three dimensions as 

oo 

= E\^){n\ (3.10) 

n 

and the phase space integral as 

/d*.^(2.)3^3fe+...+p„)-Ln/^ (3.11) 

The Wightman two point function has the spectral representation 

{Q\<P{x)<P{y)\Q) = J^^e{qo)e-^^'^-y^p{q') (3.12) 

= / dK^p(K'')A^+\x-y,K'') (3.13) 
Jo 

and the Green function is 

/•oo 

G{x-y) = / di/p{ii'^)AF{x -y]f/) 
Jo 

POO 

= ZAf{x - y; rv?) + / dp^ p{p^)Af{x - y; p^) (3.14) 
where the spectral function p{p^) in three dimensions is 

oo 

0{pW) = {27:fY.^\p-p, -P2 pr,)\{n\m\n)\' (3.15) 
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Fourier Transform equation |3.14| and make the analytic continuation to Euclidean 
space. 



= -^2+ld-'4^. (3-16) 

pZ _|_ J pZ ^ 

= m{v^) (3.17) 

The p is the multiparticle(n > 2) contribution to the spectral function. In three 
dimensions, the field strength renormalization is finite and the canonical sum rule 
holds: 

1 = y dK^p{K^) = Z + J dK^p{K^). (3.18) 

Positivity of the spectral function implies 

0<Z<1 (3.19) 

Note that the equation |3.12| is just the Fourier transform of p{q'^), so the inverse 
Fourier transform gives(choose y = 0) 

27rp(g2) = f d^xe''"'{Q\(j){x)(l){0)\Q) (3.20) 



Introduce a complete set of intermediate states p.lQ and apply the LSZ reduction 
formula to find 

27rp(g2) = GMiq)Fiq)GMiq) (3.21) 

where 

oo yn « 

^(^)=E— d<^nmi^+^M^---,Pn,q)\' (3.22) 
„ n. J 

Here we only consider the connected part of the Green function. We can use equation 
|3.21 , p.22| , and 3/7 to get the relation between a and p. Define the renormalized 
quantities as 

n + 1 

Go = Z(/^, p = Z^, Go,„^,^ = 0'(n+i) (3.23) 
and using these equations, we find 

= #^ (3.24) 



3.3 Spectral Representation of G 



From equation |3.24| , one can see that the inclusive cross section is closely related to the 
spectral function and the inverse of the Green function. The spectral representation 
of Gm{,P) in Minkowski space is 



Gm{P) 



p2 



+ 



Am? 



dp' 



— p^ + ie 



(3.25) 
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Figure 3.3: The pole and structure for g{x) 



Consider a simple complex function with the same pole and cut structure as 

g{z) = + r rfx'-^^ (3.26) 



The function has a pole at Xp and cut at a; > Xc- The pole and cut structure are 
shown in Figure |3]^. To analyze the pole structure of g~^{z), we want to know the 
zero structure of g{z). The positivity of p implies g{z) has no complex zero and no 
zero when x < Xp oi x > Xc- The only possible zero for g{^) is when Xp < ^ < Xc- 

Because g{z)g~^{z) = 1, g^^{z) should have the same cut structure as g{z). We 
can write 

c^ , r../7(x') 



g-\z) = fiz) + E ^ + / dx'^ (3.27) 

where Ci is a positive constant and 7(x) is the spectral function for g~^{x) and f{z) 
is just a polynomial function. We can simplify the calculation by considering the case 
g{z) has no zero(Ci = 0). The multiple zero case will be discussed in the Appendix. 
So the equation |3.27| becomes 

g-\z)=f{z)+ Hdx'^ (3.28) 



Because the pole of g{z) is the zero of the g ^{z), one can impose the condition 
g^^{xp) = and g{z)g'^{z) = 1 into equation ^.2(j| , ^.28l and get 



oo 



g-\z) = {z-xp)[l+ dx',-,-^!—-^ ] (3.29) 
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With the ie prescription of equation |3.25| , we can calculate the imaginary part of 



g \x) 



'^[g ^{x + ie)] = n^{x) 



^gix)] np{x) 



From equation |3.24| , we see that 7rs^cr(s) is the imaginary part of (S^jv/ . 



Hence 



Z-f{s) _ i{s) 



(3.30) 



(3.31) 



(3.32) 



To calculate the renormalization constant Z, we use g{z)g ^{z) = 1, letting the 



z ^ X. 



pi 



1 



Z J {x' — XpY 



(3.33) 



Making the analytic continuation to Euclidean space, we have the Euclidean inverse 
two point Green function 



G{p) 

Expand in power of p^, 



{p^ + rn?) 



1 



we have 



m \ 1 -\- dn 



+ pU- dK 



+ 



= ^ (m'^+p' + (■■■)/ + ■ 



(3.34) 



(3.35) 



(3.36) 
(3.37) 



In the last step, we define the m' and Z'. These quantities can be measured on 
the lattice. By taking derivatives of equations p.l6| and |3.37| , we can find relations of 
m, m' and Z, Z' as 

Z' >Z (3.38) 



m > m 
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3.4 Upper Bound of Inclusive Cross Section 



From equation |3.36| and p.37|, one can calculate yj 



dG-\p) 



Qp2 



p2=0 



1 



So we have 



4 = iW<;.= ^ 



In equation p.32| , we have defined the renormalized 7(s) as 



(3.39) 
(3.40) 

(3.41) 



From the above equation and p.40| , one can derive an upper bound on the integrated 
inclusive cross section 



ds a(s) 



Z / ds 



7(s) 



Z{ — -1 

z- 



(3.42) 
(3.43) 



Thus 1/Z' — 1, obtained from the lattice MC simulation of Euclidean G{p^), pro- 
vides a rigorous upper bound on the inclusive cross section cr(s). If 1/Z' — 1, computed 
non-perturbatively, is close to its weak coupling perturbative value, positivity implies 
that all of the multiparticle production cross sections are small. If 1/Z' — 1 is large, 
then the multiparticle production cross sections are making a large contribution. 
However, the 1/Z' — 1 only bounds the integral. Thus even if it is large, one can not 
conclude that the production cross section for any particular n is large. Note that 
from equation |3.42| and ^.38| , one can also have 



ds a{s) < ^ < 1 



(3.44) 



This implies that if the theory is strongly coupled(Z << 1), the inequality in |3.42| is 
an inefficient bound. For weak coupling case {Z' ~ Z ~ 1), the inequality in |3.42 
nearly saturated. 



IS 



Chapter 4 



Critical Energy and the Relation 
to Euclidean Momentum 

In this chapter, I will estimate the critical energy for which the multiparticle pro- 
duction turns on exponentially. First, I will calculate the multiparticle cross section 
using a statistical calculation of the multiparticle phase space integral and estimate 
the critical energy from the amplitude. Then I will compare the result with the direct 
phase space calculation for low order n. At the last, I will discuss the relation of the 
critical energy with the momentum on the lattice and determine the condition for 
the lattice simulation to be sensitive to the hypothetical large cross section for high 
multiplicity events. 

4.1 Amplitude for // ncj) 

In this section, we derive for the cross section for the production of n scalar in // an- 
nihilation and estimate the critical energy which leads to sudden exponential growth. 
Several authors 0, |, 0, H have attempted to estimate the cross section pertur- 
batively for the production of large number of 0. The most straightforward is the 
direct counting of graphs by a recursive construction, which we will outline here for 
the case of the pure 0^ without the spontaneous symmetry breaking(SSB), in which 
the recursion relation is somewhat simpler. The calculation for the case with spon- 
taneous symmetry breaking will be discussed later. For the one with SSB, we will 
follow Argyres, Kleiss, and Papadopoulo(AKP) and Voloshin's method 0, Q and do 
the calculation in 1 + 2 dimensions. 

Consider // annihilation into a single off-shell 0, which subsequently decays into 
n on-shell 0. Throughout the calculation, we only consider tree-level perturbation 
theory and neglect the contributions from // annihilation into more than just one 
off-shell (J). In the pure (p^ theory, symmetries phase, we can denote the number of 
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Figure 4.1: Diagrammatic structure of the recursion relation for 0^ interactions in 
symmetry case 

topologically distinct diagrams by a(n).(see Figure |4l|) 

, n\ ^ ^ aim) aino) airi'i) 

«H = ^^ E Sn,+n,+n,,nD{n,)D{n2)D{n,)^^^ (4.1) 

ni,n2,n2,odd TL\. Tl2- Tl^. 

The ?T,!/(?T,i!n2!n3!) counts the number of ways in which the n final particles can be 
grouped in to three sets of rii, n2 and particles, the 6 symbol accounts for that fact 
that ni + 77,2 + na = n and the factor D{ni)D{n2)D{n^) is given by the propagators 
D connecting the root vertex with the subtrees with ni,n2 and particles in each. 
Note that each a(nj) can also sub decay into another three a{nj) recursively. 

We consider the possible momentum configurations for the amplitude. We can find 
upper and lower bounds at the extreme cases of the allowed momentum configurations. 
For a generic momentum configuration, the recursion relation for the amplitude a(n) 
is too complicated too deal with. We make a simplifying assumption that leads to a 
tractable problem. We assume that for each partition of n momentum into ni+n2+n3 
momentum, the average (pi-pj) is the same for each of the momentum configuragions 
associated with rii, n2, and n^. The propagator for each internal line is the same. 



D ^{n)=q^-m^ = {pi + P2 -\ \- Pn^ - rri^ 

= np] + Y^pi ■ pj -m^ ^ j) 

= n{n — 1) {pi ■ Pj) + (n — l)m^ 

= n^{pi-pj) -n{{pi-pj) -m^) -rn^ (4.2) 

~ {n^-l){p^■p,) (4.3) 

We write 

{P^ ■ Pj) = i^m^ (4.4) 

where ^ is just a positive constant. Clearly if £^ = 1, all 3-momenta Pi must be zero 
and the phase space volume vanishes. This condition gives us a trivial lower bound 
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n = 




Figure 4.2: Diagrammatic structure of the recursion relation for mixed 0^ and 0^ 
interaction 



on the cross section. For > 1, we have 



Din, 



1 



1 



{pi + P2 + ■ ■ ■ + PiY — ni^ — 
and the recursion relation for the amplitude becomes 



(4.5) 



a{n) 



3! 



^ 5ni+n2+nz,nD{ni)D{n2)D{n-i) 



ni,n2,ng,odd 



a{ni)a{n2)a{n3) 



The solution of the above recursion equation for a{n) is 



(4.6) 



a{n) = — ^ m [n — l)n\ 



n-l 



48^2^2 



(4.7) 



Note that, we have approximated the multiparticle amplitude depending on pi, p2 ■ ■ ■ 
to an amplitude which is a function of total energy(or total number of out going par- 
ticles) only. 



For the broken symmetry case, the interaction term in Lagrange is 

A3 ,3 1 A4 4 



with 



A3 = m\/3A 



V + 



(4.^ 



(4.9) 



Use the same kinematic assumption in equation |4.3| and [4.4| , the recursion relation 
for the amplitude becomes 



n\ ^ ia[ni 

S ni,n2,n3 



tain2 



ni!(n? - 1) n2\inl - 1) n3!(n| - 1) 
ia(ni) iain2) 



2\m^ent'n:"'^"''n,\{nl-l)n2\{nl-l) 



(4.10) 
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We can make an ansatz [^: 

(4.11) 



(4.12) 

with 6(1) = 1. The /i in the case of spontaneous symmetry breaking is 

/X = = VS (4.13) 

In the special case A3 = 0, one can see that the amphtude |4.11] goes back to 
equation [4.7| . Once again, the amphtude of muhiparticle production has been ap- 
proximated by an amphtude which is a function of total energy only. 

4.2 The Phase Space 

In chapter 3, we have defined the transition probability as 

r{n, E) = ^ f d<!>n\A{p,, • • ■ ,p„)p (4.14) 
mnl J 

In last section, we have obtained an approximation of the 1 —>■ n amplitude. 

^(Pl,P2,---,Pn) >A{E) 

and the problem has been simplified to the amplitude A{E) times the volume of the 
phase space. In order to calculate the transition probability and the cross section, we 
have to calculate the volume of the multiparticle phase space in as accurate a manner 
as possible. When n becomes very large, the exact calculation becomes very difficult. 
So, we still need to find a good approximation method to calculate the phase space 
integral. One method is to use a statistical method which was developed by F. Lurcat 
and P. Mazur ||TP|] and Fermi [^. We will derive the phase-space volume based on 
their method. To simplify the calculation, we choose the center of mass frame(CM) 
for our phase space integral. 

pM = (pO,p) = (E,0) 

where 

n 



a{n) = —irn? ^'^n\{n'^ — l)h{n) 
Which leads to the following equation for the h{n) 



H») = (^5^^ ((f + v/3)" - 
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In general, the phase space integral in three dimensions has the following form: 

(4.15) 

We rewrite the phase space integral as 

jd^n = ^JdiPSUE) (4.16) 
(PSUE) = (4.17) 

and 

n „ 

Rnip) = n / d'pM - myip'r - m)SHP - Y.Pr) (4-18) 

r=l r 

We use a Laplace transformation to factorize this 

= n / d^PrC-'^P^Sipl - m2)^(p0 - m) 

= (m(^o(a))" (4.19) 

where 

m(po{a) = d^pi e-^-P'Sipf - m'^)9{p^i - m) (4.20) 

for each single particle. We can define the single particle probability distribution 
function(spd) of momentum in phase space as 

u{a,pi) = ^ 4.21 

so that 

J d^pi u{a,pi) = 1 (4.22) 

Suppose the n individual momentum pi,2- -,n are independent variables and P — 
Pi + P2 + ■ ■ ■ + Pn- Prom the statistical theory, we can define the full(all particles) 
distribution function as 

" Ud^pApl - m'')e{pl - m)e-"-^) 5\P^ - Y.rP^) 
= n ( / d^pr u{a,pr)) 5\P^ - Y.P'^) (4.25) 

We have 

. n 

y C/„(a, P) = 1, C/„(a, P) = n ^(«'Pr-) (4.26) 
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From equation |4.17| and |4.23| , we see that if we know the distribution function Un{a,p), 
we can calculate the Rn{a,p) and the phase space (PS) integral. The Un{a,p) can 
be calculated by using the central limit theorem. For large n, the sum of the n 
independent variables(P = pi + P2 + ■ ■ ■ + Pn) will have a probability distribution 
Unia,P) approaching to the Gaussian. So we can expand Un{a,P) as 



Un{a,P) = ^G{P){l + 



Ci C2 



+ 



n 



n 



where 



n 

M=0,l,2 



27r(AP) 



1 2APM J 



(4.27) 



(4.28) 



The Ci, C2 are constants and AP'^ is the variation of P^, and (P^) is the expectation 
value of P^'. In the CM frame, P^' = (P°, 0) and a ■ P = (aP°, 0). The trick is now to 
obtain a closed-form approximation for \I/(P) and then to put P'^ to its desired value 
P° = yi, P = 0. The expectation value of P^ and P^ is zero((pi) = (P^) = 0). 
Then 



AP° 
APi = AP2 



nm- 



V9o(aj 



nm 



(P2[a) 



ipi{a} 



2 \ifo{a) ) 



with the definition of y^i, v^2 as 



nr'(p2{a) 



d^Pi 5{p1 - m^)e{p'^i - m)e-'^P°p° 
d^Pi 6{pI - m')e{p^i - m)e-"*'"(p°)2 



(4.29) 

(4.30) 
(4.31) 



In the CM frame, P^ = (P'^, 0), the expectation value of P^,P^ is zero. Note that 
the asymptotic series of Un{a,P) is formally good for any P^ but useful only very 
close to the peak of the Gaussian. We therefore can tune the constant a such that 
(P°) = y/s. So, if we choose (P°) = ^/s = P°, the exponential function equal one 
and the multiparticle distribution function becomes 



Un{a,P) 



(27r)3/2(nm2)3/2 



and from equation [4.23| , we have 



RniP) 



2(mv9o(a))" e 



aE 



(27r)3/2(nm2)3/2 (^pM 



1/2 



(4.32) 



(4.33) 
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So the phase space integral becomes 



The (fo, ifi, and can be calculated directly from equation [4.20|,[4.30],and [4.31 



(4.34) 



(po{a) 
m^(y9i(a) 



Tc e 



am 
1 

— \- m] rrupo 
a 

1 A 

— + {- + my 



a 



rriLpo 



(4.35) 
(4.36) 
(4.37) 



Note that we get the equation |4.33| and [4.34| by using a statistical method which 
is accurate for large n. Although our interest is in large n, is interesting to compare 
the Rn with the exact phase space integral at n = 3. For convenience, we define l/u 
as the ratio of E and nm 

A = 1 = ^ (4.38) 

nm V lpq 

The statistical method gives us 



/7r\3/2 

U) <^-''' 

7.61 X (1 -z/) 



and the exact result is 



= -(1-^) 
2 ^ ' 

~ 4.93(1 -z/) 



(4.39) 



(4.40) 



One can see that even at small n, the form of the dependence on v is the same adn 
the coeficients are the same order of magnitude. For slightly larger n, the results 
converge rapidly. 



4.3 The Critical Energy and The relation to Lat- 
tice Momentum 

The critical energy is the energy at which the transition rate begins to grow expo- 
nentially. In section 4.1 and 4.2, we have obtained the ingrediets to approximate the 
transition rate as an amplitude squared times the phase space volume in CM frame. 

mVin.E) = |a(n)pi- / d$„ (4.41) 
n'. J 
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From equation [4.11| and [4.34| , the n particle amplitude for the broken symmetry case 
is 



and 



A{n) ~ — ^^m^(n^ — l)n\h{n) 



71-1 

2 



(27r)2-3/2 (M^ 



(/3o{a) 



(4.42) 
(4.43) 



where in broken symmetry case, A4 = 6A and A3 = ^/3. So 

b{n) = 



(i + ^r-(-i + ^r 



(4.44) 



Again, we choose the CM frame condition (P°) = P° = y/s = E. With u defined in 
4.38| and keeping just the exponential behavior, one can find the n particle transition 
rate is 

,2 



T{n,E) 



-Km 



2e 



(1 — v) X ngv{l — u) 



[n(21n(l+i/)+ln(ng)+ln(i/)+ln(l-iy))] 

2e 



where we define the coupling constant 

9 

and from E"^ ~ n^^^m^, we have 



A 



327rm 



(4.45) 
(4.46) 

(4.47) 
(4.48) 



One can see that the tree graph calculation leads to an exponentially abrupt rise of 
the inclusive cross section at some critical energy. If we fix n but send E to the 
infinity, v will vanish as well as the transition rate. If we fix E, then n can not goes 
to infinity because there is not enough energy. To find the critical energy for both E 
and n very large, we can keep the ratio v fixed at its maximum value and find the 
minimum energy of the threshold. The threshold for this blow up happens when the 
exponent is positive. So the critical energy E^, is the root of the exponent. Set the 
energy parameter e = gE/m so that 



ev = ng 



Exponential growth occurs for 



ez/2(l-z/)(l + z/)2 > 1. 



(4.49) 



(4.50) 



The function z/^(l — z/)(l + z/)^ has maximum value 0.431, so that the minimum value 
for e is e=K = 2.3198. The critical energy is 



m 

2.3198— 



(4.51) 
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In 1 + 2 dimension, the spectral integral p.l8| is convergent, so the dominant 
contributions to the spectral representation of (j (p^ ),( Equation |3 . 1 6| ) , come from p{k'^) 
with < p^- When we compute G{p^) on the lattice, the Euclidean lattice momenta 
are restricted to the Brillouin zone. 



-<V:<- (4.52) 
a a 



The restriction for will be 



Q<P^<d{^^ (4.53) 

where D is the dimension of the lattice and a is the lattice constant. Thus for the 
lattice calculation to be sensitive to cr(s) for s up to El, we must work on a lattice 
satisfying 



7r\ 



2 



Et<D[ -j (4.54) 



Dividing by m^, the requirement is 



0^>m^. (4.66) 

where rriL = ma is the mass in lattice units, which is measurable on the lattice, 
and tends to zero as the bare parameters are tunes to approach criticality(continuum 
limit). Define the correlation length in the lattice units ^ = l/m^ which should satisfy 

< — = a (4.56) 



-\/3vrm ~ 
in three dimensions. 

For a finite lattice calculation to give a good approximation to the continuum 
results in general, we should have 

a<^^ = — <^L = Na (4.57) 
m 



Use equation 4.51,4.56 one can find the condition 



TTl 

1 < 74.24^^ < ^ < (4.58) 
V 3A 



This set of conditions restricts the range of parameters of the simulation. In order 
to satisfy the critical energy threshold, the correlation length from the simulation 
should longer than ~ 42.86^. For the A/m' ~ 1 in the broken symmetry phase, this 
set of conditions require lattices of size 128^ or 256^. 



Chapter 5 
Lattice Simulation 



Monte Carlo simulations have become increasingly important as a tool in studying 
complex statistical systems and fundamental properties of quantum field theory. In 
this chapter, we will discuss how to measure Z' and m'^ on the lattice using Monte 
Carlo method. Because we are working close to the critical line, we experience the 
critical slowing down problem. I will discuss the cluster update algorithm in section 
two. We also want to know the properties of our simulation when the lattice spacing 
goes to zero. So I will discuss the continuum limit of the simulation in section three. 
In section four, we will specify the simulation method. 

5.1 Monte Carlo Simulation 

The Monte Carlo simulation on the lattice is well known. We will just follow the 
standard way [pi]] . However, it is useful to review the major features of the three- 
dimensional theory in the parameter space we use. The Euclidean action is discretized 
on an A^^ lattice with lattice constant a and periodic boundary conditions. The 
continuum mass and coupling A are related to the lattice parameters by Al = Aa and 
mi = ma. The discretized action is 



where is the lattice bare mass squared. The discretized Lagrangian can be written 
as 





[(j){n + n)(l){n) + (j){n)(l){n — n)] 



it 





n 



n 
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Since the momentum on the lattice should be inside the Brillouin zone, the discretized 
free propagator is 

Go(k) = 7 i ^ (5.3) 

{ELi^2(l-cos(aA;,))+m2} 

Note that in the simulation, we didn't distinguish different lattice momentum vectors 
which give the same momentum magnitude, i.e., \kij^k\ = \kj,k,i\ = \kk,i,j\--- etc. 

We use the Metropolis algorithm to generate a sequence of equilibrium configu- 
rations. Here we only consider the nearest neighbor interaction at each lattice site. 
The probability density is 

dP{^) = ^^Y[d<Pin) (5.4) 

n 

where the partition function is 

Z = n/ d<P{n) e"^^'^). (5.5) 



A random walk goes through the phase space with the acceptance probability 



(5.6) 



for each change of — > + (i0 at each lattice site. We measure the one-point and 
two-point functions after sweeping through the lattice. We also compute the Fourier 
transform of the field and calculate the two-point function in momentum space by 

G{k) = (0(fc)0(-A;)). (5.7) 

Then at small momentum k, the reciprocal of the two point function is 

G{k)-' = ^{m^ + P + Oik^) ■ ■ ■) (5.8) 

We then plot the reciprocal of the two-point function in momentum space against 
the reciprocal of the Fourier transform of the massless free lattice two-point function 



m = in |5.3|) . The quantities and Z' in equation 3.37 are extracted from the 



intercept and slope at the origin of this plot. To illustrate the plot, we have shown 



one example in Figure . Here we plot the Monte Carlo evaluation of the free 
propagator (A = 0) G~^{p) against the massless Gq^{p) = 4 sin^(p^)~^ = p^. 
One can see that the Monte Carlo results is the same as the known analytic result 



equation |5.3| . Note that this provids a nontrivial check of a large fraction of our 



Fortran programs. 
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Figure 5.1: Result of Monte Carlo evaluation of Gq^{p) plotted against exact analytic 
massless Gq^{p) = 4X]^sin^(p^)~^ = p ^ at 32^ lattice. Error bars on Monte Carlo 
results are too small to show up on this plot. 
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Table 5.1: Comparison of the values of the dynamic exponent z for different algorithms 
for Ising model in various dimensions [p^] . 



dimension d 


Metropolis 


Wolff 


Swendsen-Wang 


2 
3 
4 


2.167 ± 0.001 
2.02 ± 0.02 


0.25 ± 0.01 
0.33 ± 0.01 
0.25 ± 0.01 


0.25 ± 0.01 
0.54 ± 0.02 
0.86 ± 0.02 



5.2 Cluster Update Algorithm 

The Metropolis algorithm has a major drawback when working very close to the criti- 
cal point. It suffers from critical slowing down if one approaches the phase transition. 
As one approaches to the critical point, both the correlation length and correlation 



time diverge. The relation of the correlation length ^ and correlation time r are|^, ^ 

r ~ ~ (5.9) 

where the z is so called dynamic exponent and the r is measured in computer steps. 
This equation tells us how the time correlation gets longer as the correlation length 
diverges near criticality. However, in a system of finite size, the correlation length 
can never really diverge. With the lattice size L, the correlation time increase as 

T r^U (5.10) 

when close to the criticality. One can see that as we approach criticality, we need to 
have longer CPU time to get the required statistics. 

To avoid the critical slowing down, one must define appropriate nonlocal(or col- 
lective) variables and a new dynamics for driving them. Progress has been made 
with some nonlocal algorithms for both discrete spin models and continuum fields. 
Swendsen and Wang [P^ have used the Fortuin-Kastelyn percolation map for the 
Potts model to define collective coordinates that allow domains to be inverted with 
zero free energy cost. Another approach was proposed by Wolff [P^] which is essen- 
tially the same as Swendsen and Wang algorithm with some modifications. Further 
improvements are possible by introducing multigrid methods |]25[ . 

However, the algorithm for the 0^ field theory is a little more complicated. Brower 
and Tamayo||26|| proposed an algorithm for 0^ field theory based on the Swendsen and 
Wang algorithm. Here we will follow Brower and Tamayo's method but switch to the 



Wolff algorithm. Table ^]1| compare the dynamic exponent value for the Ising model 
in different dimension and algorithm. One can see that the Wolff algorithm is a little 
better than Swendsen and Wang algorithm in three dimension. 

In our simulation, the update algorithm consists basically of two parts: (1) a con- 
ventional Metropolis Monte Carlo update for the 0(x) field and (2) Wolff cluster 
identification and flipping. To identify the cluster of field, one can introduce a 
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discrete variable s^j to describe the Ising like property of the field 0, 

0(n) =s^|0(n)| (5.11) 

where the Sft = ±1. 

The detailed procedure is as follows: 

(1) . Update the 0(n) fields via a standard local Monte Carlo algorithm. In our case, 

it is the Metropolis algorithm. 

(2) . Choose a site n at random from the lattice. 

(3) . At that chosen site, n, introduce an effective Ising-type system, 

n,n+h 

and look in turn at each of the neighbor of that 0. Add the neighbor (j){n + n) 
to the cluster with the probability 

p^^^ _ ]^ _ g-|0(n)(^(n+n)|(l+SflSH+n) (5.13) 
= _ g-(|</-(n)0(n+n)|+(/-(n)0(n+n)) (5.14) 

Note that Padd = if 0^ and have opposite sign. That is, we only add the 
neighbors to the cluster when they have the same sign as the chosen one. 

(4) . From the boundary of the cluster, repeat the adding step in (3) as many times 

as necessary until there are no left in the cluster whose neighbors have not 
been considered for inclusion in the cluster. 

(5) . Flip the sign of the cluster. 

The beauty of the Swendsen-Wang and Wolff update step is that the entropy of 
the percolated cluster formation exactly cancels the surface energy to give zero free- 
energy cost. Detailed balance is guaranteed by the condition that summing out the 
percolation variables in the joint distribution, Padd-, of equation |5. 14| yields the correct 



Ising distribution of equation |5.12| . With the Wolff algorithm, our autocorrelation of 



the VEV has been reduced from ~400 to ~2. Figure |5.2| shows the differences of the 
simulation results with and without the cluster update algorithm. The simulation 
without the cluster update algorithm has six millions sweeps whereas the one with 
the cluster update only needs hundred thousand sweeps. 



5.3 Continuum Limit 



The continuum 0^ quantum field theory can be recovered from a discrete lattice theory 
by taking the limit of the lattice constant a going to zero and possibly rescaling the 
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Figure 5.2: The autocorrelation for the (0) with and without the cluster update 
algorithm. The point for the one without the cluster update algorithm is at 
about 400 (top). The autocorrelation length for the one with cluster update algorithm 
is about 2 only (bottom) . 
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field strength by a factor of Z{a). The existence of the limit of the lattice 0^ theory 
in three dimension has been proved by R. Fernandez, J. Frohlich and A. Sokal [p!l|] , 
and J, Glimm and A. Jaffe [l^. To make use of the simulation, we have to know the 
behavior of the theory as a ^ 0. 

In 0^ theory on the lattice, we have defined the = Xa and = ma which 
measured in lattice units. The correlation length in physics units is ^ = 1/m. The 
correlation length of the lattice units is = ^/a = l/mi. Consider a simulation 
on a linear lattice of A^i lattice point (in figure Ni = 8). Suppose the simulation 
with initial parameters {-Pi} gives output correlation length in lattice units (in 
figure |5.3| , E,Li =4). In the computer simulation all quantities are in lattice units; 
there is no physical lattice spacing. We are free to introduce a lattice constant ai 
in some arbitrary physical units. Then the correlation length in physical units is 
^1 = S,Li X 0.1 (= 4ai). We now want to proceed toward the continuum limit — > 0. 
We do this by a sequence of halving the physical lattice spaceing, as follows: Go to 
lattice of N2 = 2Ni[= 16 in figure p. 3D lattice points. Adjust the input parameters 



{P2} so that the output correlation length in lattice units in = 2^Li(= 8). As the 
figure shows, this may be interpreted a a2 = ai/2 and ^2 = x 02 = ^1, i-G. we 
have halved the physical lattice spacing while keeping the physical volume and the 
correlation length in physical units fixed. 



When a goes to zero, we have 



a — 







00 


rriL — 






(5.15) 
(5.16) 
(5.17) 



For the A, we have two different possibilities ||11|| : 



(1) A fixed, so that when a — > 0, Al ^ 0. This is the standard procedure in super- 

renormalizable field theory (ci < 4). From the point of view of field theory, no 
nontrivial coupling-constant renormalization is being performed; and none need 
be performed, since there is no ultraviolet divergences for the coupling constant 
renormalization in d < 4. From the point of view of statistical mechanics, the 
theory is becoming extremely weakly coupled( A^ ~ a^""^); but the effects of 
this coupling are amplified by the infrared divergences of the 0^ lattice theory 
near the critical point in c? < 4, leading to a non-Gaussian continuum limit. 

(2) Al fixed, so that when a — 0, A — 00. This is the standard procedure in the sta- 

tistical mechanical theory of critical phenomena. In this case it is for d > A that 
matters are simple: correlation function are given by a perturbation expansion 
that is free of infrared divergences. From the point of view of field theory this 
might seem surprising, since the 0^ field theory for d > 4 is perturbatively non- 
renormalizable, with horrendous ultraviolet divergence. However, it is mitigate 
by the fact that the theory is becoming extremely weakly coupled( A ~ a*^"^). 
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a J 



Mm • • • • • 

1 



a = a/2 

2 1 

Figure 5.3: The scaling of the continue hmit. When we halve the lattice constant, 
a — >■ a/2, the correlation length, will be doubled 2^^ 



On the other hand, for d < A the theory of critical phenomena is very com- 
plicated, by virtue of the infrared divergences. This is reflected in field theory 
in the fact that the theory is becoming extremely strongly coupled(A ~ a"^"^), 
causing perturbation theory to break down even though the internal momentum 
integrations are ultraviolet convergent. 

The Figure |5.4| shows the different limits in the two parameter (A, a) plane. When 
A — s> oo, that corresponds to the Ising limit. When A — 0, that corresponds to the 
Gaussian limit. When a — > 0, that corresponds to the continuum limit. The Figure 
^.5| shows the statistical mechanical approaches for < 4 and d > 4(the solid lines), 
also for the field theory approach (the dashed arrow). 

Since we want to consider a superrenormalizable theory with fixed weak coupling 
constant, we use the first alternative. Hence when a goes to zero, the A^ should also 
goes to zero. 

a — ^0, \l — ^0 (5.18) 

In the simulation, when we halve the lattice constant a a/2, we should have 
\l ^ Xl/2 and ^ ^ 2^- 

The phase diagram of the theory in terms of the parameters and A^ is sketched 
in Figure The Gaussian fixed point at (yuf^, A^) = (0, 0) is not trivial in less than 
four dimensions, as has been demonstrated by constructive field theory techniques [p^, 
[Til . The physical dimensionful coupling constant remains finite as this Gaussian fixed 
point is approached along the dashed or dotted lines sketched in Figure ^.6| . The solid 
line extending outward from the origin is the line of the second-order phase transition 
separating the broken symmetry phase from the symmetric phase. In this space of 
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Limit 



Ising Limit 



1 



Gaussian 
Limit 



Figure 5.4: The two parameters (A, a) phase diagram for different hmits 



A^oo Ising Limit 



Continuum 
Limit 




« = CXD 



Xr ,m const; d < 4 



X ,m const. 




X^,m const.; d > 4 



h=0 



Gaussian Limit 



Figure 5.5: The two parameters (A, a) phase diagram. The diagram shows the sta- 
tistical mechanical approach for both d < A and d > A. The dashed arrow shows the 
field theory approach. 
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Figure 5.6: The phase diagram for the three-dimensional, single component scalar 
0^ theory. The solid line separates the symmetric phase from the broken symmetry 
phase. The dash-dot lines represent lines of constant physics in the symmetry phase, 
and the dashed lines are lines of constant physics in the broken phase. The origin is 
the Gaussian limit and does not imply a trivial field theory for d < 4. 
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the lattice coupling constant and lattice bare mass squared, continuum physics is 
recovered by moving toward the origin along line determined by the renormalization 
group flow. Along these lines the ratio of the renormalized coupling and mass are 
constant. We want to emphasize that we are not working near the Ising fixed point, 
which in this parameterization is the point on the solid line with ^ oo. 

Other dimensionless properties {P} of the system should also remain fixed when 
we approach to the continuum limit. One expects to approach the continuum limit 
for 3> 1. If one starts far from continuum limit (^^ not ^1), the {P} will not stay 
fixed(within statistical errors) until after a number of steps have got to large enough 
^^(see figure f).3\ j. When this has been achieved, one claims that one is so close to 
the continuum limit values of the {P} that error due to lattice contamination is less 
than the statistical error of the simulation. 



5.4 The Simulation Methods 

Here we describe step by step how we carry out the simulation and analysis the data. 
Because we know the continuum limit of 1/Z' and m' exists [|T0| , p!l[] , we will try to 
halve the lattice constant a until those values settle to their continuum limits (invariant 
under a — >• a/2). The method is: 

(1) . Choose the proper values for and /ig^ as the input parameters and calculate 

the 1/Z' and m'^ by measuring the slope and intercept of the plot in section 
5.1. 

(2) . Then we check whether the correlation length satisfying the conditions in 



equation 4.58 or not 



(3) . If it satisfies the condition, then the value of 1/Z' will gives us the upper bound 

of inclusive cross section. 

(4) . If not, we halve our lattice space(go to higher energy but fix the ratio of E/n) 

and, according to the continuum limit, move to the the next point on the phase 



diagram along the dash line in Figure 5.6. In our simulation, halving the lattice 



constant and keep the volume fixed is the same as keeping the lattice constant 
fixed and double the lattice size. So we have 2N and set Ai/2 then 

find the proper value of fi^^ (actually, fi^^ /^o^/^ approximately)^ to makes 
C,L —>■ 2^L. Under this condition, we keep the dimensionless quantity X/m' the 
same. This step will guarantee we stay on the dash line in Figure ^]6| but moves 
close to the origin. 



(5). At this new scale, measure the 1/Z' and m' and check the condition [4.58| again. 



If the correlation length still can not satisfy the condition or it satisfies but the 



"'^This is because in 3-D cj)'^ theory, the bare mass square, /Iq , is hnearly UV divergent and 
propotional to uil and when a a/2, uil — > rnL/2. 
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values XjZ' and vn! do not remain the same, we then repeat step (4) until the 
correlation length can satisfy the equation [4.58| and both \jZ' and m! do not 



change anymore, (stay within estimated errors bounds) 

The last step checks the continuum limit of the simulation. If \jZ' and m' are 
insensitive to the lattice space, the value will be very close to the continuum limit (with 
some statistical errors). In the real life, we can not keep enlarging the lattice size 
because of the limit of our computing power. We can only choose the starting point 
by try and error at small lattice size iV = 32, then move up to the maximum lattice 
size N = 256. 



5.5 Verify the Ising limit of the Simulation 



From Figure |5.4| , one can see that if the A in the 0^ theory goes to infinity, it becomes 
the Ising theory. The Ising model Monte Carlo simulation have been well known in 
statistical mechanics 0. To verify our Monte Carlo simulation in three dimension 



0^ theory, one can compared the results of three dimensional 0^ at Ising limit (A 
oo) with the results from Ising model simulation. It is well known that the critical 
temperature jSc for the three dimensional Ising model is about /3c = 0.2217. The 0^ 
theory should be able to reproduce this result by letting the A ^ oo. We can rewrite 
the field theory Lagrangian |0 as 



£=-i:^#j+i:^(^f-^^) (6.19) 

where the d is the dimension. When A — oo, it should goes back to the Ising model. 
Compare with the Hamiltonian of Ising model 

n = -^Y.^ia^, (5.20) 



we see that the A — oo will force the second term in equation p.l9| to zero. So, we 
can have 

-/ii + 2d 



lim 0^ - ^11±_^ = (5.21) 

A— ►oo * A 



This suggest that 



where the ar = ±1 and 



iL-^f=±V^^-V^-r (5.22) 



P^^dl^ (5.23) 



Under this limit, the Hamiltonian of 0^ theory becomes 
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Table 5.2: The Monte Carlo simulation of 32'^ 03 theory at Ising limit(A — > oo). The 
table shows the critical temperature Pc at different A values. 



A 




Mo, 
A 


Pc 


100 


-31 


.31 


.25 


300 


-75 


.25 


.23 


1000 


-228 


.228 


.222 



= -/5E/?^J (5-24) 

which is the Ising model with the external field J = 0. 

To verify the Ising limit, we do the simulation with the A ^ 1 and calculate 
the VEV (0) for different /Iq. We start in the symmetric phase with one A value 
and change the value /Xq until the VEV is no longer zero. In our simulation, when 
A ^ 1, we get a very sharp transition from the symmetric phase to broken symmetry 
phase. So, it is very easy to identify the critical /Xq^. Table ^]2| shows the Pc from the 
simulation at 32^ with different A. The third row shows that when A = 1000, the Pc 
is very close to the value 0.2217. This give us confidence in our simulation. 



Chapter 6 

The Perturbation Calculation 



6.1 (f)'§ model and renormalization 

To compare with the simulation results, we have done the perturbation calculation 
for the VEV, {(p) / ^/m/, and 1/Z' in three dimensional 0^ theory. We use the standard 
Lagrangian 

C = ^dMx)dMx) - ^<l>lix) - ^<l>tix) (6.1) 

The field are unrenormalized canonical field, and /Iq, Aq are bare mass and coupling 
constants. To implement perturbation theory in the broken symmetry phase(/ig < 
0, {(f)) ^ 0), the field is shifted 

(0o) = Vo, 00 = Vo + 00 (6.2) 

The condition that one is perturbing about the correct vacuum is 

(0o) = (6.3) 

This condition determines V as a function of the parameters of the theory. The bare 
field and parameters are renormalized as 

0o = V^0, /io = ^m/i', V = y/ZV (6.4) 

As we mentioned before, in three dimension super renormalizable theory, there is no 
need to renormalize the coupling constant A. Hence we can set 

Ao = A, Zx = l (6.5) 

When we substitute the renormalized field in to the bare Lagrangian, and write all 
Z's as Z = 1 + SZ, one can get a renormalized Lagrangian plus terms proportional 
to V^(and containing no SZ), plus counter terms which are proportional to one or 
more 5Z's. To implement the stability condition |6]^, compute the tadpole graphs. 
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including counter terms, to a given order (number of loops) in perturbation theory, 
and expand the renormalized V 

V = C.v, (6.6) 
C. = 1 + Ci') + Cf + ■■■ (6.7) 

and adjust C?; (equivalently V) to satisfy ^]3| order by order. In the broken symmetry 
phase, we can expand the linear term in the Lagrangian(up to one loop order) 

C = -{fil + XV')V4)o (6.8) 
= -iZ^.fi'' + ZCXv'')ZC,v4> (6.9) 

-[{ZmZ - l)/i' + {Z\l - l)At;>0/ (6.10) 
+ two-loop counter terms (6-11) 



The zeroth order fixes 



= (/i^ = < 0) (6-12) 

and the one loop counter term may be written as 

£i = vSfx^^ (6.13) 
Sfx^ = {ZraZ - Z\^)Xv^ + two-loops term (6.14) 

Computation of (0) in one-loop order(Figure |6.1| ) give 

(0) = -iSXvAm + ivSfx^ = (6.15) 

where the regularized Feynman integral is 

J [Ziry p-^ — + le 

The 0^ terms in the Lagrangian gives 

A = 9090 - ^0^ + Ui^^' - '^{Z\l - 1)02 (6.17) 



are 



the one-loop counter terms(for — zS) generated by equation |6.17 

CT = i{Z - l)q^ + i6fi^ - i{Z'^Cv - 1)^^ (6.18) 

The one-loop self-energy diagrams are shown in Fig. 11. The tadpole diagram plus 
the tadpole counter term add to zero by our previous choice of Sfi"^. Hence, the 
renormalized self-energy function is 

- iS(g2) = -ilSX'^v^Immiq^) + iSZq^ - i{Z^Cv - (6-19) 





(X) 



Figure 6.1: The one-loop tadpole diagram and the counter term. 



o 



Figure 6.2: The one-loop self energy diagrams. 
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To determine the renormalization constants, we can impose the on-shell renormal- 
ization conditions 

d_ 



E(g2 = m') = 0, ^E(g2^ 



= 



(6.20) 



This give us the on-sheU renormahzation constants and S(y 

9A 



6Z 

sCv 



9A 
47rm 



-- + iln\/3 
3 2 



(6.21) 
(6.22) 
(6.23) 



Alternatively, the Z' and m' are defined by the zero momentum renormalization 
conditions 

d_ 



= 0) = 0, ^E(5^) 



= 



(6.24) 



which give 



5Z' = 

K = 

5Z' = 



3A 



327rm 
21A 



327rm 

39A 3A 

+ 



327rm 27rm^ 



(6.25) 
(6.26) 
(6.27) 



6.2 One-Loop calculation for 



We calculate the vacuum expectation value (VEV) to one-loop order and compare to 
the simulation results. In the simulation, we measure the bare VEV and 1/Z' and 
ml. In three dimensions, (0o) has dimension -^/mass. Hence, the quantity {(po)/\/m/ 
is dimensionless and we can look for its continuum limit. In the broken symmetry 
phase, the VEV will is 



So we have 



(00 



(0o) 



m 



V — 



2A 



(6.28) 



Im' / 9A 13\ 
2A I ^4^48/ 



(6.29) 
(6.30) 
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6.3 The calculation for ^ to 0{X'^) 

In chapter 3 (and Appendix A), we have derived a bound on the inclusive cross section. 

ds a{s) <^-l (6.31) 

The lattice Monte Carlo simulation p[rovides a nonperturbative result for the right 
hand side of equation |6.31| . To see if this implies any large nonperturbativecontribu- 



tion, we have to compare it to a perturbative case of the equality. (The C introduced 
in appendix A is zero in the weak coupling limit.) 

^ - l) = / ^(^) = ^ = E -^^""^ (6-32) 



The index (A^) indicate the order of A^. So that 

Z' Z 



\ - 1 = - = A'(i) - bZX^^^ + ;f + . . . (6.33) 



where 

=Y.\ds af ) (6.34) 
Here we have to sum over all final states with n outgoing particles. 

C^n{s) = / d^rX^mi+M. ■■■,Pn, P)\' (6.35) 

where the phase space integral is 

l^^- = h.l ^^(^^ - • ■ ■/ ^Avl - m')i2nmP - Y^p.) (6.36) 



The lowest order, 0(A), contribution has only one diagram. Figure shows the 
lowest order tree diagrams. So that 



1 1 

2 

9m^X 



d^2 (6.37) 



If we choose the CM frame, the two particle phase space integral will be 



Hence the contribution of order 0(A) is 



/<i4, = ji= (6.38) 
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Figure 6.3: The one-loop Feynman diagrams for n = 2(0*1 




aa 



bb 



+ 




+ 




+ 




Figure 6.4: The tree diagrams for n = 3(^2)- 



By using equation 3.21 , the second term in equation 6.33 will be(in order of 0(A )) 

2 



5^(l);t'« = - 



3A 



327rm 



;-8 + 61n3) 



(6.40) 



The Feynman diagrams and the counter terms for the third term in equation 
|6.33| are shown in Figure |6.3| and Figure |6.4[ Figure |6.3| shows the one-loop Feynman 



diagram with n = 2{G—d>mm/l) and Figure |6^ shows the tree level Feynman diagram 

3 

for n = 3(0*4). The cross section of order is 



1 1 



27rs2 2 
1 1 

27rs2 3! 



<l>323?[0'|0'|] 



(6.41) 
(6.42) 
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We can write 



The one-loop Feynman integrals (with r 

.9y2At / 6 



in Figure are 



Us) 



6TT 



X2 



In 



his) 

his) 
his) 
C.T. 



[y2 

yi 

.9V2X 



-1 + X2)x2l/2 / 



In 



-2 + X2 -2^2^ 

xly2 



(-1 + xi)xiyi J 



2m 



8tt 
.9V2X' 



In 



y/s + 2rri 
^ — 2m ^ 



X 21n3 



6TT 



9V2At 
t— X 21n3 = /, 



in 



.9V2X^ 



in 



4 



where xi, X2, and i/i are 



Xi 
X2 

yi 



2 + 2^/f^ 

4 - r 



2-2^/7^ 
4 - r 



~l~ 1 



(6.43) 



(6.44) 

(6.45) 

(6.46) 
(6.47) 
(6.48) 



(6.49) 

(6.50) 
(6.51) 



The two particle phase space integral is the same as before(equation |6.38| ). Finally, 
after calculating the s integral, one gets 



(2) 



4m^ 



ds af^ 
3A 



inm 



X 0.496728 



(6.52) 
(6.53) 



The Figure |6]^ shows the n = 3 case of 0(A ). 



(2) 



and 



'^3 



1 1 



27rs2 3! 



/ ds 

rf$3 |0'4|' 



(6.54) 

(6.55) 
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From Figure |6]J, the four point function is 

(pi = Taa + Tbb 



the Taa is just the tree diagram in Figure |6.4| aa. 

Taa = -iQX 



The tree level diagram in Figure 6Abb has six permutations. 



T 



bb 



1 1 
+ 



1 



Si2 — m"^ S23 — m? S31 — m? 



(6.56) 
(6.57) 

(6.58) 
(6.59) 



where 

^ij — -\- kj) 

is the sum of the momentum of two of the final particles. We make the same simplifi- 
cation that was made on the calculations of the multipicle phase space ie. we replace 
each ki ■ kj by the average (fcj ■ kj). Then 



(si2 - m^) = {s23 - rr?) = (sgi - m^) = - 
Hence, the four point function square can be simplified to 



|^^(s)|2 = 36a2 1 + 54 



— + 729 — 

s s 



(6.60) 



(6.61) 



Since the four point function does no longer depend on the individual momentum of 
each final particle, the phase space integral becomes easy. We can use the recursion 
relation [4.17| in chapter four to calculate the space space integral. We find 



/ 



167r 



So the contribution from becomes 



(2) 



ds 
3A 



1 1 



3m 



3.62) 



2713^6 

-I X 0.201276 



(6.63) 
(6.64) 



Combining equation |6.39y6.40|j6.53| , and |6.64| gives 



1 3A 1 

— = IH X - 

Z' Svrm 4 



3A 



X 0.206944 



(6.65) 



The mass m in above equation is physical mass under the on-shell scheme. How- 
ever, we only measure m' from Monte Carlo simulation. In order to compare with 



48 



the Monte Carlo results, one should use m! in the perturbation series. The relation 
between m! and m can be calculated from two point function. 



Go(p^) = Z'G'{p') = ZG{p'), 



and 



,/2 



Go^(p^ = 0) = ^ = iK + S(p^ = 0)) 



So up to one-loop order, we have 



Z' 



ZZI = + =0)) 



l + 6Z'-5Z + —^iO) 



1 + A 



From equation |6.7CI| , we have 



m 



A 



Airm V 24 2 



Plug into equation |6.65|, we get 



3A 1 

X 



3A 



?7rm' 4 V 87rm' 



X 0.201214 



(6.66) 
(6.67) 

(6.68) 

(6.69) 
(6.70) 

(6.71) 
(6.72) 



(6.73) 



Chapter 7 

Analysis and Results 



In this chapter, I will discuss the results from the Monte Carlo simulation and compare 
the results with the perturbation results. 



7.1 The Results From the Simulation 

The results for the broken symmetry phase are given in Table |7]1|. In going down 
from line to line in the Table we have N 2N, C,l 2^^ (approximately). 
This corresponds to holding the physical lattice size L = Na fixed while halving the 
physical lattice spacing a. Here we keep the ratio L/C,l ~ 5.6. This procedure will 
keep on a line of "constant physics" as approach to the continuum limit (Figure |5.6|) . 
Note that the dimensionless quantities A/m', ^ — 1, and {(f)o) / ^/m/ keep the same 
values, within errors, while we are halving the lattice space. This shows that these 
values are insensitive to the lattice constant a and approaching to the continuum 



limit. The Figure |0, IT^, and 0| plot the A, 4t and ^4= in Table \fA\. The plot 



shows the insensitivity of lattice size for both those quantities. The dashed lines are 
the fitted lines which show small slopes of the plot. However, it is within the error 
bar. 



Table |7.2|J7.3| , and |7.4| shows the finite size effect of the simulation. We have 
2N, ^l( approximately) while the input parameters /Xq^, are held fixed. 

This corresponds to doubling the physical length of the lattice while keeping the lattice 
a fixed. This check the sensitivity of the output results to lattice size corrections. The 
Tables show small discrepancies at smallest lattice size. However, we found that if we 
keep the ratio N/C,l > 3, the finite size correction can be ignored. In Table |7.1| , we 
keep the ratio N/C,l ~ 5.6 so that we can ignore the finite size effect while we are 
moving along the RG trajectory. 
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Table 7.1: The results of our simulations on different size lattices(A^). The output 
correlation length = l/m'^^, the dimensionless ratios X/m', and 1/Z' are given. 





Xl 


Ar3 




m' 


1 

Z' 


\frnJ 


-0.161 


0.2 


32=^ 


5.7(1) 


1.13(2) 


1.008(2) 


0.78(1) 


-0.078 


0.1 


643 


11.3(4) 


1.13(4) 


1.012(3) 


0.80(2) 


-0.0384 


0.05 


1283 


22.11(1.0) 


1.11(5) 


1.012(5) 


0.80(2) 


-0.01905 


0.025 


2563 


45.4(3.0) 


1.13(8) 


1.010(8) 


0.79(3) 



Table 7.2: The finite size effect with fixed input parameters /Xq^ = 0.078, = 0.1. 





Xl 


iV3 


iL 


\_ 

m' 


1 

z> 




-0.078 


0.1 


323 


12.9(4) 


1.29(4) 


1.0025(10) 


0.73(1) 


-0.078 


0.1 


643 


11.3(4) 


1.13(4) 


1.012(4) 


0.80(2) 


-0.078 


0.1 


1283 


11.2(4) 


1.12(4) 


1.020(5) 


0.80(2) 


-0.078 


0.1 


2563 


10.8(8) 


1.08(8) 


0.98(2) 


0.78(3) 



Table 7.3: The finite size effect with fixed input parameters = 0.0384, Xl = 0.05. 





Xl 


ATS 




m' 


1 

Z' 


(■^0) 


-0.0384 


0.05 


323 


21.1(1.6) 


1.06(8) 


1.0002(12) 


0.71(1) 


-0.0384 


0.05 


643 


22.1(1.0) 


1.10(5) 


1.003(2) 


0.77(2) 


-0.0384 


0.05 


1283 


22.1(1.0) 


1.10(5) 


1.012(5) 


0.80(2) 


-0.0384 


0.05 


2563 


21.9(1.1) 


1.09(6) 


1.015(9) 


0.79(3) 



Table 7.4: The finite size effect with fixed input parameters yUg^ = 0.01905, Xl = 
0.025. 





Xl 


iV3 




m' 


1 

Z' 


iM 

\frn! 


-0.01905 


0.025 


643 


41.5(6.5) 


1.04(16) 


1.002(3) 


0.78(2) 


-0.01905 


0.025 


1283 


45.1(3.8) 


1.127(97) 


1.0011(35) 


0.77(2) 


-0.01905 


0.025 


2563 


45.4(3.0) 


1.13(8) 


1.010(8) 


0.79(3) 
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Continuum Limit 
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Figure 7.1: Continuum limit for the The plot shows the insensitivity of lattice 
size for 

m 
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Figure 7.2: Continuum limit for the J7. The plot shows the insensitivity of lattice 
size for 
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Figure 7.3: Continuum limit for the The plot shows the insensitivity of lattice 
size for 
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7.2 Discussion and the Conclusion 



From Table [7.1| , on the 256^ lattice, we have 



^, = 45,4 > Ii^ip = 37.9, (7.i; 



This value satisfies condition |3.43| . From the table, the l/Z' — 1 is 

^ - 1 = 0.01 ±0.008 (7.2) 

At 95% confedence level, the upper bound of the inclusive cross section is 

y ds o{s) < 0.026 (7.3) 

One can see that the numerical result shows no evidence for large multiparticle pro- 
duction at high energy in weakly coupled 0| theory. To gain further insight into the 
significance of this results, we compare it to the perturbative calculation. 

1 3A 1 / 3A \^ , , 

— = 1 + X - - X 0.201214 7.4 

Z' Sirm' 4 {Snm' J ^ ' 

From Table |7. 1| , at the = 256^ lattice, we have A = 1.13. So the perturbation 
result is 

^ = 1 + 0.1349 X - - (0.1349)2 x 0.20121 (7.5) 

= 1 + 0.03372- 0.00366 (7.6) 
= 1.03006 (7.7) 



and the 1/Z' in Table |7]T| shows 

^ = 1.010 ±0.008 (7.8) 

Zj 

we can see that there is little room for multiparticle contributions to the inclusive 
cross section. 



We can also do a perturbative calculation of (^o) /ym' and compare to the lattice 
Monte Carlo result. The result is striking! The perturbation shows 

> i^(l + J^}f\ (7.9) 



For ^ = 1.13, the result is 



m' V 2A V 47rm' 48 



°^ - 0.665(1 ±0.219) = 0.811 (7.10) 
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The VEV from our simulation(A^ = 256) is 
0.022. So the simulation results is 



]) = 0.118 and from Table ^ 



m 



0.795 ± 0.097 



(7.11) 



In conclusion, our results in equation |7]^ rule out the possibility of large nonper- 
turbative effects coming from high multiplicity production amplitudes in the weakly 
coupled 0| theory. The comparison with perturbation results from above leaves little 
room for even small nonperturbative effects. 



Appendix A 

Multiple zeros case for the 
function g~^{z) 



In chapter 3, equation p.271 leaves open the possibihty of zero of g{z){po\e of g (z) 
We deal explicitly with the case of one zero of g~^{z). For the multiple zeros case, 
one only need to replace the single C term by J^Ci- So the equation |3.27| becomes 



g-\z) = fiz) + 



C 



+ 



dx' 



x' — z 



(A.l) 



Figure |A.lj shows the possible function of g{z) with one zero, note taht C is positive. 



g ^{xp 



Because the pole of g{z) is the zero of g ^(2;), one can impose the condition 
and g{z)g~'^{z) = 1 into equation p.26| , |A.l| and get 



g \z) = {z-Xp) il 



+ 



dx' 



7(x') 



{x' — z) {x' — Xp) ^ 



(A.2) 



In multiparticle production at high energy, x » Xc, equation [A.2| has a pole on 
the cut and will give us the results the same as equation |3.30| and p. 32 . 



To calculate the renormalization constant Z, we use g{z)g ^{z) = 1 and take the 



1 
Z 



1 + 



c 



dx 



{x' - XpY 



+ 



(A.3) 



After analytic continuation to Euclidean space, the inverse two point Green function 
becomes 



Expand in power of 



1 + 



+ / dK^ 



1 



2 4 

— + — + 

^2 ' ^4 ' 



(A.4) 



(A.5) 
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Figure A.l: The function of g{z) with one zero. 



we have 



G-\p) = m^il + 



C 



+ / dK^ 



C 



+ 1 + — + / rfK 



^4 ' / ^4 



+ 



7 (m'^+p^ + (■••)/ + 



(A.6) 
(A.7) 



From equations |A.6| and [A.7|, we have 



1-1 + -^+ Id.'^-^ 



(A.8) 



From equations p.41| and |A.8| , the upper bound of the integrated inclusive cross 
section becomes 



ds ais) = Z ds 



\z' e^ 



(A.9) 
(A.IO) 
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< {^-l) (A.11) 

Because the C is a positive number, the bound remains. It is just a less effecient 
bound. 
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